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Abstract. As a part of the ISOPHOT CIRB (Cosmic Infrared 
Background Radiation) project we have searched for point-like 
sources in eight fields mapped at two or three wavelengths be- 
tween 90/im and 180/itn. Most of the 55 sources detected are 
suspected to be extragalactic and cannot be associated with pre- 
viously known objects. It is probable, also from the far-infrared 
(FIR) spectral energy distributions, that dust-enshrouded, dis- 
tant galaxies form a significant fraction of the sources. 

We present a tentative list of new extragalactic FIR-sources 
and discuss the uncertainties involved in the process of extract- 
ing point sources from the ISOPHOT maps. Based on the an- 
alyzed data we estimate the number density of extragalactic 
sources at wavelengths 90/im, 150/mi and 180/im and at flux 
density levels down to lOOmJy to be l«10 5 sr -1 , 2-10 5 sr~ 1 , 
and 3 TO 5 sr _1 , respectively. 

Strong galaxy evolution models are in best agreement with 
our results, although the number of detections exceeds most 
model predictions. No-evolution models can be rejected at a 
high confidence level. 

Comparison with COBE results indicates that the detected 
sources correspond to <;20% of the extragalactic background 
light at 90/im. At longer wavelengths the corresponding frac- 
tion is ~10%. 

Key words: Galaxies: evolution - Galaxies: starburst - Cos- 
mology: observations - Infrared: galaxies 



1. Introduction 

The cosmic infrared background (CIRB) consists in the far- 
infrared of the integrated light of all galaxies along the line 
of sight plus any contributions by intergalactic gas and dust, 
photon-photon interactions (7-ray vs. CMB) and by hypothet- 
ical decaying relic particles. A large fraction of the energy 
released in the universe since the recombination epoch is ex- 
pected to be contained in the CIRB. An important aspect is the 
balance between the UV-optical and the infrared backgrounds: 
what is lost by dust obscuration will re-appear through dust 



emission in the CIRB. Some central, but still largely open, as- 
trophysical problems to be addressed through CIRB measure- 
ments include the formation and early evolution of galaxies, 
and the star formation history of the universe. 

The primary goal of the ISOPHOT CIRB project is the de- 
termination of the flux level of the FIR CIRB. The other goals 
are the measurement of its spatial fluctuations and the detec- 
tion of the bright end of FIR point source population contribut- 
ing to the CIRB. The full analysis of the data from the DIRBE 
(H auser et al. |1998| ; Schlegel et al. |l998[ ) and FIRAS (Fixsen et 
al. 1998) experiments indicated a CIRB at a surprisingly high 
level of ~ 1 MJy sr" 1 between 100 and 240 /im. Preliminary re- 



sults had been obtained already by Puget et al. (1996). Lagache 



et al. (1999) detected a component of Galactic dust emission 
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associated with warm ionized medium and the removal of this 
component lead to a CIRB level of 0.7 MJy sr _1 at 140/im. 

Because of the great importance of the FIR CIRB for cos- 
mology these results definitely require confirmation by in- 
dependent measurements. ISOPHOT observation technique is 
different from COBE: (1) with relatively small f.o.v. ISOPHOT 
is capable of looking at the darkest spots between the cirrus 
clouds; (2) ISOPHOT has good sensitivity in the important FIR 
window at 120 - 200 /im; (3) with the good spatial and spectral 
sampling ISOPHOT gives the possibility of recognizing and 
eliminating the emission of galactic cirrus. 

In the ISOPHOT CIRB project we have mapped four low- 
cirrus regions at high galactic latitude at the wavelengths of 90, 
150, and 180 /im (see Fig. [[]). Through this multi-wavelength 
mapping we will try to separate the cirrus component and con- 
firm the detection of sources at neighboring wavelengths. In 
addition, we have performed absolute photometry in several 
filters between 3.6 - 200 /im at the darkest spots of the fields. 
This photometry will be used ( 1 ) to secure the zero point for the 
maps at 90, 150, and 180 /im, and (2) to determine the contri- 
bution by the zodiacal emission using measurements of its SED 
at mid-IR wavelengths where it dominates the sky brightness. 

This paper presents the first step in the analysis of the 
ISOPHOT CIRB observations. Here we will concentrate on 
the data reduction and the study of the point sources (galax- 
ies) found in the FIR maps. The source counts determined in 
the FIR are important for the study of the star formation his- 
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tory of the universe and for the testing of the current models of 
galaxy evolution. 

With recent observations at infrared and sub-mm wave- 
lengths it has become obvious that star formation efficiencies 
de rived f rom optical an d UV observ ations only (e .g. Madau et 
al. |1996[ Steid el et a l. |1996| , |l999fc Cowie et al. |1996| , |l997 



1998; Hu et al. 1998) underestimate the true star formation ac- 



tivity at high redshifts because the correction for dust extinction 
is unknown (e.g. Heckman et al. 1998). 

IRAS has shown that in the local universe about one third of 
the luminosity is emitted at infrared wavelengths. In starburst 
galaxies the fraction can be much higher as most of the starlight 
is absorbed by dust and re-radiated in the infrared. In extreme 
sources like the hyperluminous galaxy IRAS 10214+4724 the 
energy spectrum peaks around 100/im in the rest frame and 
more than 90% of the energy is emitted in the infrared and sub- 
mm regions. The emission maximum moves further towards 
sub-mm with increasing redshift, causing optical studies to se- 
riously underestimate the true star formation activity. If the dust 
content is high enough the objects can remain completely un- 
detected at optical wavelengths. 

With ISO and new sub-mm instruments like the SCUBA 
bolometer array (Holland et al. 1999) it has become possible 



to study the star formation history of the universe at infrared to 
sub-mm wavelengths (for reviews see e.g. Hughes et al. 1998, 



1999). Due to the negative A'-correction the observed flux den- 



sities will not depend strongly on the redshift and it is possible 
to detect more distant galaxies (e.g. van der Werf 1999; Guider- 
doni et al. |1997|). 



1999; Lilly et al. 1999 Blain et al. 1999a) have shown that star 



Recent studies (e.g. Dunlop et al. 1994; Omont et al. 1996 
Hugh es et al. |1997|, |1998| ; Stiavell i et al |l999t Abraham et al. 



formation activity remains high at z > 1 . 

In observations with t he SC UBA instrum ent at 450/im and 
850m n (e.g. Hugh es et al . |1998| ; Barge ret al . |1998| ; Smailetal. 
1997| ; Blain et al. |l999<j ; Eales et al. |1999| ; Lilly et al. |1999| ) 
galaxies have been detected up to redshifts z— 5. Compared 
with galaxies seen in optical surveys the objects have higher 
dust content and the star formation rates are an order of mag- 
nitude higher. The surface density of the detected sources ex- 
ceeds predictions of no-evolution models by at least one order 

, Eales et al. 



of magnitude (Smail et al. 1997 



1999;Barger et al. 



1999). The number of sources detected by Eales et al. (1999) 
3 mJy accounts for -20% of the CIRB de- 
1998). Similar results were 



at 850/im above 
tected by FIRAS (Fixsen et al 



obtained by Barger et al. (1999). At the level of 0.5 mJy the 
sources contain most of the sub-mm CIRB ( Smail et al 



1999; Blain et al. 1999b) 



Kawara et al. 



observerved the Lockman Hole at 



95/im and 175/im using ISOPHOT. The number of sources 
found was at least three times higher than predicted by no- 
evolution models. The conclusions of the FIRBACK (Puget et 
al. |1999| ) and ELAIS (Oliver et al. |2000| ) projects are similar 
and at 175/im sources with S'i/>120mJy account for —10% of 
the CIRB detected by FIRAS. 



In this article we study the number density of extragalac- 
tic sources and their contribution to the FIR background radia- 
tion using observations made with ISOPHOT. The data consist 
of maps made at wavelengths 90/im, 150/im and 180/im, and 
for some smaller areas at 120/im. The total area is close to 1.5 
square degrees. Most of the regions have been observed at three 
wavelengths (90/im, 150/im, 180/im) some at two wavelengths 
(120/im and 180/im). Both the galactic foreground cirrus emis- 
sion and the emission from typical extragalactic objects will 
reach their maxima within or near the observed wavelength 
range. In particular, we will be able to determine the cirrus 
spectrum for each region separately. 

We have developed a point source extraction routine based 
on the fitting of the detector footprint to spatial data. The 
method is different from those used in most previous studies 
where the source detection algorithms have concentrated on the 
variations (off-on-off) of the detector signal as function of time. 
Our analysis will therefore be independent of and complemen- 
tary to the previous results. 

2. Observations and data processing 

The observations were p erformed with the ISOPHO T pho tome- 
ter (Lemke at el. |1996|) aboard ISO (Kessler et al. |1996|). The 



maps were made in the PHT22 staring raster map mode using 
filters C_90, C_120, C_135, and C_180 with reference wave- 
lengths at 90/im, 120/im, 150/im and 180/im, respectively. In 
this paper we study eight fields that cover a total of 1 .5 square 
degrees (see Fig. [j]). Two fields were observed at two wave- 
lengths while the rest have observations at three wavelengths, 
the details are shown in Table [l]. The CI 00 detector used for 
the 90/im observations consists of 3 x 3 pixels each with the 
size of 43.5"x43.5" on the sky. The C200 detector used in 
the rest of the observations has a raster of 2x2 detector pix- 
els, 89.4" x 89.4" each. 

We have selected regions with low surface brightness. 
Some maps have redundancy i.e. the observed pixel rasters 
partly overlap each other. In four maps observed with filter 
C_90 the raster step is larger than the size of the detector, lead- 
ing to incomplete sampling. 

The data were first processed with PIA (PHT Interactive 
Analysis) program versions 7. 1 and 7.2. Special care was taken 
to remove glitches caused by cosmic rays since these might 
be erroneously classified as point sources during later analysis. 
The flux density calibration was made using the FCS (Fine Cal- 
ibration Source) measurements (FCS1) performed before and 
after each map. Generally the accuracy of the absolute calibra- 



19g3 tion is expected to be better than 30% (Klaas et al. |1998[). The 



calibration was normally applied to observations using linear 
interpolation between the two FCS measurements. 

The data reduction from the ERD (Edited Raw Data; de- 
tector read outs in Volts) to SCP (Signal per Chopper Plateau; 
signal at each sky position in units V s^ 1 ) was performed also 
using the so-called pairwise method (Stickel, private comm.). 
Instead of making linear fits to the ramps consisting of the de- 
tector read-outs one examines the distribution of the differences 
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between consecutive read-outs. The mode of the distribution is 



estimated with myriad technique (Kalluri & Arce 1998) and 



Simulated measurements were also used to study the effects 
of imperfect fiat fielding (Appendix B.4). These do not cause 



is used as the final signal for each sky position. This process- 
ing was done in batch mode. Compared with the previous PIA 
analysis there were some calibration differences which were 
possibly due to the different drift handling of the FCS mea- 
surements. For this reason the final surface brightness values 
in the new AAP (Astrophysical Applications Data) files were 
rescaled using the results of the previous interactive analysis. 
The subsequent analysis was carried out with both data sets but 
no significant differences were found in the results. The pair- 
wise method is, however, believed to be more robust against 
glitches and in the following the results are based on the data 
reduced with this method. 



3. The point source detection 

The point source detection was performed in two steps using 
data processed to the AAP level with PIA and the pairwise 
method (see Sect. |J). The data consists of surface brightness 
values with error estimates and were flat fielded using special 
routines (see Appendix B.l). 

In the first phase each surface brightness value was com- 
pared with the mean and the standard deviation estimated from 
other measurements within a fixed radius. The radius was set 
to about three times the size of the detector pixel. Values ris- 
ing above the local mean surface brightness by more than 0.7 a 
were considered as potential point sources. 

In the second phase a model consisting of point source and 
a constant background was fitted into each region surrounding 
the candidate positions. In the fit the footprint matrices were 
used to calculate the contribution of the point source to the ob- 
served surface brightness values. The free parameters of the fit 
were the source flux density, the two coordinates of the source 
position, and the background surface brightness. Details of the 
fitting procedure are given in Appendix B.2 



3.1. Completeness and false detections 

As a result of the fitting procedure we obtained estimates for 
the flux density, coordinates of the potential point source and 
a value for the local background brightness. The fitting routine 
provides formal error estimates for all free parameters. In ad- 
dition, we have calculated the standard deviation of the surface 
brightness values inside the selected region. Assuming that the 
flux density values are normally distributed we can calculate 
for each source candidate the probability, P, that the detection 
is not caused by background noise. 

The completeness of the point source detection and the 
probability of false detections were studied with simulations 



(see Appendix B.3 ). The number of false sources was found to 



somewhat exceed the expected number of l-P per map pixel 
and this was taken into account in setting the selection criteria. 
A confidence limit of P=99% (corresponding to ^2.3 a) was 
used to discard uncertain detections. 



significant errors in the source counts. Cosmic ray glitches and 
short term detector drifts may still lead to spurious detections 
primarily on the CI 00 detector. Glitches have been removed 
during the standard data reduction and any remaining anoma- 
lies should be reflected in the error estimates calculated for the 
surface brightness values. Large error estimates will reduce the 
number of false detections (see ^|). 



For the purpose of source counts (see Sect. |5.1[ ) we will use 
an additional criterion based on the ratio between the source 
flux density, S v , given in Jy and the background rms noise, 
(Tb g , given in units of Jy per pixel, 



P : 



S v 



(1) 



The parameter p is not directly related to the probability ob- 
tained from the footprint fit and can be used as safeguard 
against false detections. In Sect. 5.1 we will use a limit 
p > po =10.5. In the maps fluctuations are typically below 
O.lMJysr -1 which, in the case of the C200 detector, corre- 
sponds to <Tbg ^0.02 Jy per pixel. Our p-criterion implies thus a 
typical detection limit of 200 mJy. Because of the smaller pixel 
size of the CI 00 detector, —44" instead of ~89", the detection 
limit is lower by a factor of four. Note that since the source 
flux is always distributed over several map pixels p cannot be 
interpreted directly as a a limit. 

The value p n =10.5 was chosen based on simulations: at 
this limit the detection rate is ^70% for sources which fulfill 
the previous criterion of P >99%. The number of false detec- 
tions is less than one per 200 map pixels. For the C200 detector 
the expected number of false detections within the mapped ar- 
eas is SilO i.e. less than the probable number of true sources 
above the p limit that remain undetected. Because of the larger 
number of map pixels in the 90pm maps there can be more false 
detections than undetected true sources. However, the relative 
error in the source counts should remain well below 50%. 

3.2. Multiwavelength confirmation of source detections 

Since the source counts might be slightly overestimated even 
after applying the additional p-criterion constraint we have 
constructed another source list based on multi-wavelength de- 
tections. A source is accepted only if there are detections at two 
wavelengths, each at 99% confidence level (see above), and the 
spatial distance between these detections is <80". The p cri- 
terion is not used. If a map contains several possible sources 
within that distance the one with the highest estimated proba- 
bility is selected. With the number density of detections at the 
99% confidence level we can estimate the probability of finding 
by accident a source within 80" radius of any given position: 
5% for the 90^m and 150/im maps and below 10% even for 
180/im. 

We can calculate a rough estimate for the number of false 
detections as follows. For any individual detection we have an 
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Table 1. The eight fields studied here. The columns are: (1) the name of the field used in this paper, (2)-(3) coordinates of the 
centre of each field, (4)-(5) galactic coordinates of the field, (6) area of the map, (7) the number of raster positions observed over 
the map, (8) the step between adjacent raster positions in the staring mode mapping and (9)-(12) the integration time for each 
filter. The distance of adjacent scans was in all cases identical to the raster step used along the scan line 
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5x3 


135 




36 
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1% probability that it is not caused by a real source. Empiri- 
cally, the probability of finding another detection at a different 
wavelength within 80" radius is less than 10%. This increases 
the confidence in the first detection from 99% to 99.9%. In the 
source counts for 150/im and 180/im the total number of spu- 
rious detections can therefore be estimated to be ~1. Because 
of the larger number of pixels in the 90/im maps the number 
of false 90/im sources could be up to ^6. However, the false 
detections are made only close to the background noise level 
where the number of undetected true sources can easily ex- 
ceed this. Counting only detections with multiwavelength con- 
firmation will underestimate the true number of sources. The 
probability of finding two spurious detections within 80" is 
0.001% and the confidence level of a source detected at two 
wavelengths (instead of a detection at one wavelength only) 
equals a 4.2a detection. 

Examples of such multiwavelength detections are shown in 
Fig. H The positions of four sources in the field EBL26 are 
overlaid on the 90/im, 150/im and 180/im maps. Sources 5 and 
6 are detected only at 150 and 180 /mi (see also Table 5. 1 ). The 
leftmost source is a 90/im detection only and is therefore not 
in Table 5.1. It is probably caused by a glitch as can be seen 
from the SRD data (Signal per Ramp Data i.e. signals, V s _1 , 
derived from individual detector integration ramps) in Fig. ||a. 
The signals are shown for four pixels at five positions (given 
as a function of time) centered on the position closest to the 
fitted source position. Source 7 was detected only at 150/im and 
180/im while at 90/im it remained below our detection limit. 
The SRD data for this source are also shown in Fig. ||. One must 
remember that the detection procedure was not based on SRD 
data and in Fig. || only part of the relevant data are shown. SRD 
signals were, however, used to visually estimate the possible 
effect of glitches on the source detections (see^). 

For most of the area of the 15x15 position 180/im map 
NGPN no observations at other wavelengths were available. 
In order to confirm the detections the data were divided at the 
ERD level into two parts. The first data set contained the first 
ramps of each measurement and the second data set the rest of 



the ramps. These data sets were used to create two maps and 
sources were detected and accepted according to procedures 
described above. The maps are not completely independent. 
A large glitch that affects several ramps may have influenced 
both data sets. The risk for false detections is therefore larger 
than in other cases. Some of the detections could, however, be 
confirmed with observations of the 32x4 position NGPN map 
which partly overlaps the square 15x15 position 180/im map 
(see Table |l|). 

4. Discrimination against cirrus knots 

We will check the possibility that some of the sources detected 
are small scale cirrus structures (cirrus knots). 

4.1. Cirrus confusion noise 

We estimate the contribution of cirrus to the sky confusion. The 
total noise, ototah in the source flux densities is estimated from 
the standard deviation of the flux densities. This is obtained by 
fitting a point source to the position of each measurement (each 
pixel and raster position) while keeping the source positions 
fixed i.e. effectively convolving the map with the detector foot- 
print. The total noise consists of measurement noise, a meas ., 
from the detector, and the sky confusion, cr s ky, caused by real 
sky brightness variations, 



-'total 



'sky 



(2) 



The measurement noise, in flux density units, is derived from 
the error estimates obtained for the surface brightness values 
from the SRD data. The surface brightness map was modi- 
fied by adding corresponding amount of normally distributed 
noise. The fitting procedure was then repeated. The rms dif- 
ference between the flux densities obtained from the two maps 
gives an estimate of measurement errors on the flux density 
scale. Table || summarizes the results for the fields EBL26 and 
NGPS. The sky confusion is obtained (using Eq.||) from the to- 
tal and measurement noise estimates. The measurement noise 
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Fig. 1. The positions of the observed FIR fields overlaid on the all-sky map composed of the DIRBE 60, 100, 140 and 240/im 
Annual Average Sky Maps. The 180/im ISOPHOT maps are shown for four regions. Besides dark areas the fields contain, on 
purpose, also areas with faint cirrus 



was seen to be roughly equal to the sky confusion. In order to 
estimate the contribution of cirrus to the sky noise we have used 
the approach of Gautier et al (1992). 



actly the same as in our detection procedure but should give 
accurate estimates of the expected cirrus confusion. 



Herbstmeier et al. (1998) estimated for the large NGPN 
Gautier et al. (1992) calculated the confusion noise due to 180/im map a fluctuation power P =2.3 x 10 3 Jy 2 sr _1 at the 

scale of 4'. The fluctuations are mostly due to cirrus emission. 
With the dependence P ~ d 2 derived by Herbstmeier et al. 
we obtain a fluctuation amplitude of 18Jysr~ 05 at d =1.5' 
which corresponds to the size of the detector beam. According 



infrared cirrus for different observation strategies. We use these 
results to estimate cirrus contamination in the case of circular 
aperture with diameter d, immediately surrounded by a refer- 
ence annulus of the same width. This configuration is not ex- 
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Fig. 3. SRD signals for two sources shown in Fig. ^ a source detected at 90/im only (a) and the Source number 7 (b). The 
signals are shown for four pixels as function of time for five consecutive raster positions. For each pixel the plot is centered on 
the position closest to the average positions of the detections made at different wavelengths. The signals from this position are 
between the dashed lines. The pixel numbers are shown in the figure. The first source was detected at 90 /im only while the 
second one was detected at 150/^m and at 180/im 



to the tables of Gautier et al. ( 1992 ) the estimated flux density 
fluctuations due to cirrus emission are 8.0 mJy. In this field all 
detected sources have flux densities exceeding lOOmJy. It is 
therefore very unlikely that they could be caused by cirrus. 

The surface brightness attributed to cirrus is obtained 
by subtracting the CIRB given by Fixsen et al. (1998), 
0.82 MJy sr -1 at 180/im, and the zodiacal light according to 



Leinert et al. (1998). In the field, EBL26, which has the largest 



surface brightness among our fields, the average cirrus sur- 
face brightness is ^1.4MJysr _1 compared to ^0.8 MJy sr -1 
in the NGPN field. With the relation P ~ < B > 3 (Gautier et 



al. 1992 ) between the fluctuation power and the mean surface 
brightness we can estimate that in EBL26 the cirrus fluctuation 



level is ^43mJy. All sources are above ^150mJy, i.e. above 
the 3.5 er level. 



The results of Herbstmeier et al. (1998) show that at the 
scale of the C100 beam size, d ~ 45", the expected cirrus fluc- 
tuation amplitude is clearly below 10 Jy sr~ 5 for all our 90/im 



maps. According to Gautier et al. ( 1992 ) this corresponds to a 
flux density of 4 mJy which is again clearly below the flux den- 
sities of even the faintest sources at 90 /mi in the present study. 

In conclusion, it is clear that cirrus is not the dominant fac- 
tor in the sky confusion noise. The flux density values derived 
for cirrus contamination are small compared with the faintest 
sources detected here. Because of the non-gaussian nature of 
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MJy/sr 
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15.51 



14.71 



14.10 



13.68 



MJy/sr 

9.33 



8.41 



7.66 



7.13 



6.80 



MJy/sr 

6.53 



6.14 
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5.54 



5.34 



Table 2. The estimated total rms noise, Ototab m tne fields 
EBL26 and NGPS compared with the estimated measurement 



"20 s l h 18 m l h 17 n 
RA (2000.0) 

Fig. 2. The positions of four detections in EBL26 overlaid on 
the 90/im, 150/im and 180/im maps. The leftmost source was 
detected only at 90/im while the other sources (numbers 5, 6 



and 7) were detected at 150/im and 180/im. See Table 5.1 for 
parameters of the numbered sources 



the cirrus fluctuations (Gautier et al. 1992) it is, however, pos- 
sible that our source lists contain some cirrus knots. 

We have, however, also the advantage of having observa- 
tions at different wavelengths which offers an additional means 



noise, cr m cas., and the sky confusion, 
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36 


36 
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41 


40 



for the discrimination against spurious detections caused by 
galactic cirrus. 

4.2. Cirrus spectra 

Using observations made at different wavelengths we can de- 
termine the cirrus spectrum in each region. For each 180/im 
measurement the corresponding surface brightness values at 
the other wavelengths were calculated using weighting with 
a gaussian with the approximate size of the C200 detector 
beam. Linear fits were performed to surface brightness values 
at 90/im or 150/im plotted against the 180/im values and the 
slopes were used to derive the cirrus spectrum (see Juvela et al. 
2000). 

|7l The cirrus 



The spectra obtained are shown in Figs, 
spectrum can be determined most reliably in regions with clear 
surface brightness variations which is the case for fields EBL26 
and NGPS. It is also, nevertheless, possible to determine the 
spectra for all other fields. Bright individual sources were al- 
ways removed from the data but in the VCS and VCN fields 
the results may be influenced by the emission from the nearby 
galaxy, NGC 5907. 

4.3. Source spectra vs. cirrus spectra 

The cirrus spectra were compared with the spectra of the de- 
tected sources. The relationship between the cirrus surface 
brightness spectrum and the "source" spectrum due to a cir- 
rus knot depends on the size of the cirrus knot: (1) If the cirrus 
knot is small compared with the CI 00 detector pixel, the knot 
spectrum is the same as the surface brightness spectrum. (2) If 
the knot size is >45" some part of the 90/im flux is included 
into the background, and the source flux density at 90/im drops. 
(3) If the cirrus knot is large compared even with the pixel size 
of the C200 detector also at the longer wavelengths only part 
of the total flux density contained in the cirrus knot will be de- 
tected. In that case there will be probably no detection with the 
CI 00 detector. 

Simulations showed that for gaussian cirrus knots with 
FWHM ~90" we detect with the CI 00 detector less than half 
of the total flux density. For cirrus knots with FWHM^180" 
the ratios F(90/im)/F(150/zm) and F (90/rni)/F(180//m) drop 
to Si0.25 of the cirrus surface brightness values. 
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If the source spectrum is fiat compared to the cirrus 
surface brightness spectrum we can conclude the source is 
not a cirrus knot. Otherwise, this alternative cannot be ex- 
cluded with certainty. However, when the flux density ratio 
F(90^m)/F(150/xm) or F(90/xm)/F(180/mi) is clearly larger 
than one fourth of the corresponding ratio in the cirrus surface 
brightness spectrum it is improbable that the source is due to 
a cirrus knot because this would require that the knot is much 
smaller than 180". From the slope of the cirrus power spec- 
trum as functi on of spatial freque ncy, (G autier et al. 1992 ; Low 
& Cutri 1994 ; Herbstmeier et al. 1998 ) and the very low level 
of cirrus emission in most of the fields, we estimate that it is 
very improbable that there are many small (^90") cirr us knots 
strong enough to be detected as sources (see Sect. |4.l| ). 

In Figs. |[]7] we show for our four main fields the pre- 
dicted spectra caused by cirrus knots together with those ac- 
tual sources that were detected at 90/im and at some longer 
wavelength. The cirrus spectra have been drawn for two cases 
assuming either that the cirrus knot is smaller than the CI 00 
detector pixel (solid line) or that it is of the same size as the 
beam of the C200 detector pixels (dashed line). 

The source spectra are almost without exception flat com- 
pared with the steeper of the two cirrus spectra and the sources 
are not likely to be caused by cirrus. Furthermore, in many 
cases the relative 90/im flux is higher than what is possible even 
for very small cirrus knots. No sources were rejected from our 
sample based on their SED shape. 

5. Source counts and identifications 

5.7. Cumulative source counts 



Table 5.1 contains sources detected at two or three wavelengths 
with an individual detection probability P >99% at one wave- 
length (see Sect. 3.2). The coordinates are averages, weighted 
with the estimated confidence of detections at different wave- 
lengths. The flux densities and their error estimates are ob- 
tained from the fitting routine. If the field was mapped at three 
wavelengths but the source was detected only at two wave- 
lengths, then we quote an upper limit based on the background 
surface brightness variations, S — 10.5<7bg where <7b g is given 
in units of Jy per pixel. The detection probability correspond- 
ing to this limit is 97%. 

The surface density of sources is estimated by dividing the 
number of sources counted at each flux density level by the to- 
tal effective map area. With "effective" map area we mean that 
portion of the total map area from which such sources could 
have been detected. For example, in the field VCN faint sources 
cannot be detected close to the bright galaxy NGC5907. The 
maps around Mrk314 and ZWIIB were originally observed in 
order to study a known object and since these bright sources 
are not the result of random selection they are excluded from 
the source counts. 

As a first approximation we have calculated the effective 
map areas according to the faintest detected source. At each 
flux density level the area was taken to be the sum of the ar- 
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Fig. 4. Average source spectrum of the fields VCN and VCS 
(filled squares) and the cirrus spectrum determined in the same 
fields (open squares). The shaded area indicates the possible 
range for cirrus spectra. At 90/im the highest cirrus flux is ob- 
tained when the knot is small compared with the CI 00 pixel 
size and the lowest flux density corresponds to a knot with size 
similar to the C200 detector pixels. The upper cirrus spectrum 
and the average source spectrum have been fitted with modi- 
fied Planck Law, v 2 B v , with temperatures shown in the figure 
(upper frame). Spectra of sources in fields VCN and VCS de- 
tected at 90/im and at one or two longer wavelengths. The aver- 
age source spectrum has been calculated without including the 
source with the highest 90/im flux density (lower frame) 



eas of those maps where sources with equal or lower flux den- 
sities were detected. The procedure is likely to underestimate 
the surface density of faint sources since it ignores the back- 
ground surface brightness variations within the maps. The de- 
tection of just one source with a flux density S v increases the 
corresponding area A{S U ) by the the area of the whole map. 
On the other hand, for statistical reasons the faintest detected 
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Table 3. (See next page) Sources detected at two or more wavelengths. The coordinates are weighted averages of the positions 
fitted at different wavelengths. For sources undetected at one wavelength we give upper limits which, according to simulations, 
correspond to the flux density of a source with detection probability 97%. Quality flags, q, are based on visual inspection of 
the SRD data and the scale is from (obvious glitch) to 4 (clear source). Counts c are the results from our detection procedure 
based on the SRD data (see^). Asterisks mark cases where a detection at one wavelength could be associated with two different 
detections at another wavelenth. The confidence level of each detection at one wavelength is 99% while the confidence level for 
a source with two detections at different wavelengths is better than 99.99% 



=1 
O 
CO 



in 



10" 




10 - 



^source 22. 6K 



T — 15 QK 

1 cirrus 



I .11 




t I I I I I I L 



100 150 
A [/iiri] 



200 



o 

CO 



GO 



10" 



>, 10" 

!-3 



10 e 



T =24 OK 

1 source ^-^iv 




^cirrus 20. 2K 



III 




j i i i i i i i i i_ 



100 150 
A [jiirri] 



200 



Fig. 5. Cirrus spectrum and source spectra in the field EBL22 Fig. 6. Cirrus spectrum and source spectra in the field EBL26 
(symbols as in Fig. |5] (symbols as in Fig. |] 



source can sometimes be brighter than the faintest theoretically 
detectable source but this should not lead to significant under- 
estimation of the areas at any flux density levels. 

For 90/im maps observed with a 180" raster step there are 
gaps between the rasters. We have chosen not to correct for this 
effect but note that it may lead to an underestimate of the source 
densities by some tens of per cents i.e. the effect is comparable 
with the calibration uncertainties. 



The cumulative source densities obtained at 90/^m, 150/im 
and 180/im are shown as histograms in Fig. ^|. Two sets of 
sources were used in deriving these curves. The first set con- 
sists of all detections (i.e. P >99%) and no confirmation was 
required at a different wavelength (dotted line). The second set 
contains only those sources that were confirmed by a detection 
at another wavelength (solid line). 

The first set may contain a number of false detections. The 
second set gives a conservative estimate for the true number of 
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0.10(0.04) 


2/2 




5 


02 28 19.7 


-25 56 55 


0.10(0.03) 
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<0.14 




0.14(0.04) 


3/8 


EBL26 


1* 


01 16 38.9 


+02 30 45 


0.36(0.06) 
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Fig. 7. Cirrus spectrum and source spectra combined from the 
fields NGPN and NGPS (symbols as in Fig. f§ The plotted av- 
erage source spectrum has been calculated without including 
the source with the highest 90/im flux 



sources. Since the areas used in deriving the surface densities 
depend also on the selection criteria applied, the sample with 
the larger number of sources does not necessarily lead to higher 
source density. The results obtained from the two sets are very 
similar for 150 and 180/im. At 90/im the counts based on de- 
tections confirmed at 150/im or 180/im are, however, below 
the other estimates. Most of the difference can be explained by 
statistical uncertainties. There are only three confirmed 90/im 
sources brighter than 400 mJy and at lower flux levels where 
the number of sources is larger the estimates are in fair agree- 
ment with each other. The change in the ratio between the two 
estimates as the function of source flux is not statistically sig- 
nificant. Some of the difference may be caused also by the 
source properties. Sources detected at 150/im are likely to be 
seen also at 180/im (and vice versa) while more of the 90/im 
sources remain unconfirmed at the longer wavelengths. 
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Fig. 8. Cumulative source counts at 90/im, 150 fim and 
1 80/im. The solid and the dotted lines represent values obtained 
by estimating the areas according to the faintest sources de- 
tected in a map. The dotted line gives all sources detected and 
the solid line sources detected at more than one wavelength. 
The dashed curves represent cases in which both the the source 
selection and the area determination were based on the local 
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We have derived a third set of cumulative source densi- 
ties by selecting sources based on the ratio p between the flux 
density and the background surface brightness variations (see 



Sect. 3.1). For each flux density level all sources with p > po 



were selected and no confirmation at other wavelengths was 
required. The area corresponding to a given flux density level 
was obtained by first calculating the local standard deviation 
of the surface brightness values around each observed position 
and then integrating the total area with noise below 1/po times 
the source flux density. 

The limiting value, po =10.5, was selected based on simu- 
lations (see Sect. 3.1) which indicate that close to the detection 



limit the number of false detections is still clearly less than the 
expected number of true sources. In the measurements the typ- 
ical background noise is 0.02 Jy per pixel and the given p value 
corresponds to source flux density limits 45 mJy and 200 mJy 
for CI 00 and C200, respectively. For higher flux densities the 
number of false detections drops rapidly while for lower flux 
densities more of the true sources are either rejected based on 
the p criterion or are not detected at all. At the quoted flux lev- 
els the cumulative source counts are higher than the number of 
the false detections which will therefore not affect the results 
significantly. However, especially at 90pm, the number of false 
detections can exceed the number of undetected true sources 
and the source counts could be slightly overestimated. 

The cumulative source densities obtained with this third 
method are drawn in Fig [| with dashed lines. These results are 
based on the integrated area of the regions where background 
noise is low enough for a source with given flux density to be 
detectable. Since the area determination and the source detec- 
tion are based on similar criteria the errors caused by an in- 
correct p-limit tend to cancel out. At the bright end the results 
agree with the earlier histograms since no sources are rejected 
and the corresponding areas converge towards the total area 
mapped. The differences are more pronounced at faint flux den- 
sities. When the area corresponding to a faint flux density limit 
was determined based on the faintest observed source the area 
was likely to be overestimated and the source densities under- 
estimated. When the area was determined by the local proper- 
ties of the background brightness and the sources were selected 
using a related criterion the results should be more reliable. On 
the other hand, the applied p limit excludes many of the fainter 
(but more uncertain) sources that were included in the previous 
counts and the curves cannot be extended reliably to equally 
low flux densities. The values obtained for 150pm and 180pm 
below 150 mJy are probably only indicative. 

5.2. Cumulative flux densities 

Fig. ^ shows the cumulative flux densities, F V (S V ), i.e. the sur- 
face brightnesses due to all sources brighter than a given flux 
density S v . Compared with source counts the flux density val- 
ues are more sensitive to bright sources. 

Systematic calibration errors would cause the curves to 
be shifted horizontally. Our simulations show that the source 
flux densities depend also on the detection process. The fitted 



source positions can be slightly displaced from the true posi- 
tions towards a direction where the background noise produces 
the highest values. Therefore it is possible that the faintest 
flux densities are overestimated. The errors are, however, only 
noticeable close to the detection limit and should never be 
more than 10%. Furthermore, several maps were observed with 
partly overlapping rasters. The improved sampling should lead 
to more reliable flux densities. The largest uncertainty is there- 
fore in the absolute calibration itself. 

5.3. Association with known sources 

Table | lists IRAS Point Source Catalog (PSC) and Faint 
Source Catalog (FSC) sources within the mapped regions. 
ISOPHOT sources detected within a -1' radius of the IRAS 
positions are also shown. 

In the regions studied there are four sources with IRAS de- 
tections at 100pm. In EBL26 the PSC source at l h 17 m 10.6 s 
+2°17'10" has been detected also with ISOPHOT. The flux 
density obtained at 90pm, 2.6 Jy, is somewhat higher than the 
IRAS 100pm value of 2.2 Jy. The FSC source at l h 18 m 05.0 s 
+1°58'59" with 0.84 Jy flux density at 100pm has not been de- 
tected and the quoted 90pm upper limit is 1 lOmJy. The reason 
is that the 90pm map was incompletely sampled and the IRAS 
position lies between the observed rasters. 

In VCS the IRAS source at 15 h 15 m 53.3 s +56°19'47" is 
the bright galaxy NGC5907. The source is extended and the 
ISOPHOT map is too narrow for the estimation of the total flux 
density. The ISOPHOT values in Table [i] are obtained by fitting 
of a point source within an area with 5' radius. The flux den- 
sity at 90pm is slightly lower than the IRAS value. In the case 
of Mrk314 the flux densities at 120pm are comparable to the 
IRAS 100pm values. There are additionally a number of IRAS 
sources with upper limits at 100pm. In all cases the ISOPHOT 
upper limits derived at 90pm are much lower although in the 
case of maps observed with 180" raster steps the upper limits 
are not always accurate. 

In Table ^ we list objects from the Simbad database that are 
within the mapped regions. The list contains all sources iden- 
tified as galaxies or quasars. Other types of objects (e.g. radio 
sources) are listed only if they are close to a FIR detection. 

6. Discussion 

6.1. Comparison with other ISOPHOT source counts 

At the flux level of lOOmJy the following source densities are 
obtained (dotted line in Fig. ||): 1.4xl0 5 sr~\ 2.5xl0 5 sr _1 
and S.SxlO^r- 1 at 90 pm, 150 pm and 180 pm. In Figjl^ 
these results are compared with results from other ISOPHOT 
projects. 



At 90pm the results of Kawara et al. (1998) are some 30% 



higher than our source counts while the results of the ELAIS 
survey (Oliver et al. 2000| ) and Linden- V0rnle et al. (2000) are 
lower than our counts. The ELAIS counts based on observa- 
tions of 1 1.6 square d egrees at 90pm, are a factor of three be- 
low the Kawara et al. (1998) value. 
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Table 4. IRAS Faint Source Catalog and Point Source Catalog sources within a ~1' radius of the positions listed in Table 5.1 



Table 5. Sources from the Simbad database together with FIR detections within ~1.5' radius. The table contains all galaxies and 
QSOs inside the regions mapped and all other objects close to sources detected in the ISOPHOT maps 



The calibration adopted in the ELAIS project is based on 
DIRBE and it was found that the PIA analysis resulted in higher 
surface brightnesses (Efstathiou et al. 2000). This is consistent 
with our findings. The difference between the DIRBE surface 
brightness values and our data calibrated with PIA version 7.3 
is ^30% (Q). With DIRBE calibration our source count points 
at 90/im move towards smaller flux densities (see Fig. [k]) and 
they would be in good agreement with the ELAIS results. 

At 180/im we can compare our source counts with Kawara 
et al. ( |1998[ ) and Puget et al. ( |1999| ) observations at 175/tm. 
Our counts are almost two times higher than the Kawara et al. 
results but compatible with Puget et al. (1999). 



180/im are 0.13 MJy sr -1 , 0.67 MJy sr" 1 and 0.82 MJy sr" 1 , 
respectively. Adopting these values our sources would con- 
tribute already a significant fraction of the CIRB, especially 
at 90/im where ^20% of the CIRB could be attributed to the 
detected sources. For sources brighter than 100 mJy the corre- 
sponding fraction at 150 /im is slightly below and at 180/im 
slightly above 10%. Recalibrating the 140/im and 240/im 
DIRBE observations using the results of the cross-calibration 



between FIRAS and DIRBE (Fixsen et al. 1998) would, nat 



urally, move the D IRBE points closer to the FIRAS values 
(Fixsen et al. 1998 ). As discussed in |a|, our calibration appears 
to be closer to the FIRAS than the DIRBE scale. 



6.2. Comparison with galaxy models 

At 150/tm and 180/im the counts are much higher than pre- 
dicted by no-evolution models (e.g. Guiderdoni et al. 1998); at 
180 /im the difference is a factor of five (see Fig. |l0|) 



6.4. Comparison with galaxy spectra 



In the evolutionary model E by Guiderdoni et al. (1998) 
both the star formation rate and the relative number of ULIRGs 
increase with z. The model has been found to be in good agree- 
ment with extragalactic background light measurements in both 
optical and FIR. Our source counts exceed, however, the model 
predictions at all three wavelengths (see Fig. |l0| ). 

Franceschini et al. (1998) have presented similar mod- 



els which include contributions from two galaxy populations: 
dust-enshrouded formation of early-type galaxies and late-type 
galaxies with enhanced star-formation at lower redshifts. The 
predicted source counts at 1 80/im are higher than in the Guider- 
doni et al. model E and the model is therefore in better agree- 
ment with our results. 

6.3. Extragalactic background radiation 

Measurements of the CIRB in the wavelength range of the 
present observations have been recently published based on 
DIRBE and FIRAS observations of the COBE satellite. After 
the removal of interplanetary and galactic foreground sources 



(Kelsall et al. |1998| Arendt et al. |1998[ ) the level of CIRB 
detected by DIRBE was found to be 25±7nWm" 2 sr _1 at 
140 /i m and 14±3nWm -2 sr -1 at 240/im (Hauser et al. 
1998). These numbers correspond to 1.1 MJy sr -1 . According 



to Fig. |9J the sources found in this study represent some 5% of 
the CIRB, the exact number depending on which source counts 
are applied. 

Fixsen et al. ( 1998| ) have reported results based on three 
different methods used to subtract the galactic foreground 
emission. From the analytical representation of the average 
spectrum the surface brightness values at 90 /im, 150 /im and 



In Fig. 1 1 we compare the average of the source spectra pre- 
sented in Figs. |]-|7] with the spectra of the galaxies Arp 193 and 
NGC4418. In the sample of luminous infrared galaxies pre- 
sented by Lisenfeld et al. (2000) Arp 193 has the lowest and 
NGC4418 the highest estimated dust temperature. 

In the rest frame the SEDs of luminous infrared galaxies 



reach maxima between 60/tm and 100/im (Silva et al. 1998 
Devriendt et al. 2000; Lisenfeld et al. 2000| ). The spectra of 
our FIR sources are relatively flat in the observed wavelength 
range and peak typically above 90/im. This is consistent with 
most sources being at redshifts 0.5 SSz^l. 

The emission maximum of norm al spiral galaxies is also 
close to 100/im (e.g. Silva et al. 1998 ). However, in the case of 
a spiral galaxy a FIR detection at the level of 0. 1 Jy would cor- 
respond to a visual magnitude brighter than 16 and the optical 
counterpart should be visible. The lack of visual counterparts 
indicates that most of our sources are likely to be distant lumi- 
nous infrared galaxies. 

We have conducted a follow-up study in the VCN region. 
A ^2x2' field surrounding the positions of the FIR sources 



VCN 1 and 2 in Table 5.1 has been observed in V- and /-bands 



(images provided by J.C. Cuillandre, CFHT) and in K— and 
J-bands (images taken by P. Vaisanen). 

Our photometry reveals a number of red sources with 
I — K il 3. These are potential candidates for being lumi- 
nous infrared galaxies (LIRG) and could be the counterparts 
of our FIR detections. In Fig. [l2| we show the flux densities of 
these sources together with the ISO detections. In the same fig- 
ure the spectral energy distributions of two luminous infrared 
galaxies, Arp 220 and Mrk231, are drawn as they would be 



seen at redshift 1.0 (Ivison et al. 1998). The figure shows that, 



in principle, any one of the red NIR sources could be responsi- 
ble for our ISO FIR detections. An unambiguous identification 



14 



M. Juvela et al.: Far Infrared Extragalactic Background 



t- 


i 


i iii 

90yu,m 


1 


A 








" i\ 














- 


.. v . 










v.": 


■ — ...... 

v. '- - - 




i i 




" I 

i iii 


U. ... i- 

L "l 

i 



10 



S y (mJy) 



10' 



150/xm 




10 



S l/ (mJy) 



10" 




10 



S^mJy) 



1(T 



Fig. 9. Cumulative flux densities, F v , at 90/i, 150 fim and 
180/im in units of surface brightness. The curves correspond 
to samples obtained with the same methods as in Fig. || 
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Fig. 10. Comparison with other ISOPHOT counts and models 
of galaxy evolution. Our source counts at lOOmJy an d 200 mJy 
are shown together with th e results of Puget et al. ( 1999 ; cir- 
cles), Kawara et al. ( |l998| ; tr iangle s), Oliver et al. ( |2000| ; di- 
amonds) and Linden- V0rnle (2000; stars). We present our re- 
sults as horizontal lines that indicate the difference between 
the DIRBE calibration and the adopted ISOPHOT calibration. 
At 90/im the DIRBE calibration results in lower flux densities 
and at 150/im and 180/im in higher flux densities. The pre- 
dictions of model E of Guiderdoni et al. 1998 are shown with 



solid lines and the evolutionary model of Franceschini et al. 



1998 with a dashed line. Dotted lines show predictions of no- 



evolution models (90/im: Guiderdoni et al. 1998; 150/im and 
180/im: Franceschini et al. 1998) 



is not possible. Normal elliptical or spiral galaxies, as shown in 
the figure, are excluded. 
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Fig. 11. The average of the source spectra shown in Figs. 
a nd the two-temperature model spectra of Lisenfeld et 
al. ( |2000 i for the luminous infrared galaxies Arp 193 and 
NGC4418. Dashed lines show the spectra of the two galaxies 
shifted to £=1.0. The flux density scale is arbitrary 

7. Conclusions 

We have searched for FIR point sources in raster maps ob- 
served with the ISOPHOT CI 00 and C200 detectors at wave- 
lengths between 90/im and 180/im. The total area covered is 
~1.5 square degrees. Most of the FIR sources detected are pre- 
sumably IR galaxies which, due to the negative if-correction, 
can be observed at redshifts z > 1. A comparison of the SEDs 
of sources detected at 90, 150, and 180/im with cirrus spec- 
tra shows that for most sources an explanation in terms of cir- 
rus knots can be excluded. Based on the number counts of the 
sources we can conclude: 

- We have found 55 FIR sources that, due to the multi- 
wavelength confirmation, correspond to detections with 
high confidence level (^4cr) 

- We have derived a FIR source density of ~60 sources per 
□° at lOOmJy level 

- The source density is much higher than predicted by no- 
evolution galaxy models; at 180/im the excess is close to a 
factor of five 

- The source counts are in agreement with models where the 
star formation rate and the relative number of ULIGs in- 
creases strongly with z, e.g. the counts are slig htly h igher 
than predicted by model E of Guiderdoni et al. ( 1998 ) 

- At 150/im and 180/im the combined flux of detected 
sources accounts for ~10% of the CIRB intensity as de- 
rived from the COBE FIRAS data; at 90/im the fraction is 
over 20% 
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Fig. 12. Flux densities of sources with I-K > 3.0 (filled 
squares) together with the flux densities of the two ISO detec- 
tions (at 90/im only an upper limit) in a field within the VCN 
region. The spectra of the brightest sources with arbitrary I- 
K within 50" of the FIR source positions are shown as open 
squares. The spectral energy distributions of two luminous in- 
frared galaxies, Arp 220 (solid line) and Mrk231 (dashed line), 
normalized to the 150/im observations, are drawn for z ~1.0 
(Ivison et al. 1998 ). For comparison, template spectra of an el- 
liptical galaxy (dotted line) and of the spiral galaxy M5 1 (dash- 
dotted) are shown for z = (Silva et al. Il998b 
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Appendix A: Calibration comparisons 

The results presented in this paper are based on the calibration 
performed with the onboard FCS. This has been estimated to be 
better than 30% for extended sources brighter than 4MJy sr _1 



(Klaas et al. 1998). Because of the very low surface bright- 
ness and the correspondingly low FCS power, it is interesting 
to compare the ISO calibration scale with the corresponding 
DIRBE values. 

For the comparison of absolute surface brightnesses we 
took the DIRBE ZSMA (Zodi-Subtracted Mission Average 
Maps) data around each ISO field. Zodiacal light was added 
to the DIRBE values according to the model given by Lein- 
ert et al. ( |1998| ) using the A — Aq and the ecliptic latitude of 
the ISO observations. We could also have used the Weekly Av- 
eraged Sky Map data observed during those weeks when the 
solar elongation is similar as in the ISO observations. Because 
of the somewhat lower noise we chose to use the the ZSMA 



16 



M. Juvela et al.: Far Infrared Extragalactic Background 



Table A.l. Comparison of the DIRBE and ISO flux density 
scales for surface brightnesses. The columns are: (1) the field 
name, (2) wavelength of ISO observations, (3) mean surface 
brightness of the ISO map, (4) the ratio of the DIRBE and 
ISO absolute surface brightnesses. For column 4 the numbers 
given in parentheses give the dispersions in the DIRBE surface 
brightness values interpolated to the wavelength of the ISO ob- 
servations (see text). Note that for EBL22 we have chosen to 
correct the calibration according to existing absolute photome- 
try measurements 



Field 


A 


< S > 


Sdirbe/Siso 




Om) 


(MJy sr" 1 ) 




(1) 


(2) 


(3) 


(4) 


EBL26 


180 


6.0 


1.44(0.30) 




150 


7.9 


1.40(0.29) 




90 


15.8 


0.79(0.05) 


NGPS 


180 


3.1 


1.29(0.17) 




150 


3.7 


1.24(0.15) 




90 


5.7 


0.65(0.06) 


NGPN 


180 


1.8 


1.67(0.07) 




180 1 


2.4 


1.28(0.03) 




150 


2.5 


1.43(0.07) 




90 


5.5 


0.57(0.02) 


NGPN & NGPS 


180 


2.2 


1.44(0.23) 




150 


3.3 


1.31(0.18) 




90 


5.6 


0.62(0.06) 


EBL22 


180 


2.3 


2.25(0.09) [1.21(0.05)] 








1.73(0.07) 




150 


2.0 


2.00(0.07) [1.24(0.05)] 








1.69(0.06) 




90 


0.8 


0.80(0.02) 








0.85 


ZWII 


180 


12.6 


1.28(0.03) 




120 


9.6 


1.59(0.04) 



1 the larger 180[im map 

I ] using default responsivities 

italic entries: corrected according to absolute photometry measure- 
ments 



data. Both DIRBE and ISO were colour corr ected assuming a 
v 2 B u sp ectrum and T dust =18 K (Arendt et al. |1998[ Schlegel et 
al 



1998 ). DIRBE data were interpolated to the wavelength of 
the ISO observations. The average ISO flux density weighted 
with the DIRBE beam was compared with the DIRBE value. 



The surface brightness ratios obtained are shown in Table A. 1 
in column 4. The main uncertainties are due to the large DIRBE 
beam and the large noise of the DIRBE 140 /im and 240 /im ob- 
servations in these faint surface brightness regions. The given 
error estimates correspond to the total dispersion in the surface 
brightness values. These estimates are probably more realistic 
than those obtained by combining the error estimates calculated 
for the mean surface brightness values. 

At 90/xm the DIRBE values are -30% lower than the ISO 
values. On the other hand, there is an indication that the ISO 
C200 values are -30-40% lower than the DIRBE ones. 



At 140/im and at 240/im the relative error estimates given 
for ZSMA surface brightness values exceed in many cases 
50%. When the temperature was fixed to 18 K the interpolated 
values at 150 /im and 180 /im are mostly determined by the se- 
lected temperature and the 100 /im DIRBE values which have 
significantly smaller error estimates. An increase of Tdust by 
1 K increases Sdirbe/Siso at 90 /im by —1% and decreases 
it at 150/im and 180/im by about 10%. On the other hand, if 
the exponent in v a is decreased from 2.0 to 1.0 the discrep- 
ancy between ISO and DIRBE scales would increase at 180/im 
by —20%. If all sources of uncertainty are taken into account 
our results do not indicate a difference between the ISO and 
DIRBE flux density scales exceeding —30%. 

In the case of EBL22 there was an unusually large differ- 
ence, by a factor of —0.6, between the values obtained with 
default responsivities and those calibrated with the FCS mea- 
surements. The fact that the FCS heating power was clearly 
below the range for which calibration tables exist may have 
contributed to the large discrepancy. Comparison with DIRBE 
supports the higher surface brightness values obtained with the 



default responsivities (see Table A.l). Absolute photometry at 
one position in the field gives values closer to the default cal- 
ibration. Since the FCS calibration is more reliable in the case 
of the absolute photometry we have re-scaled the maps to agree 
with the absolute photometry. For all the other maps the differ- 
ences between the FCS and the default responsivities remained 
below 30%. 



Appendix B: Details of the point source detection 

B.l. Flat fielding 

Normally the flat fielding was done by calculating the ratios 
between each detector pixel and the average of other measure- 
ments within a small area around it. Linear fits provided the flat 
fielding correction factors as function of the surface brightness. 

In some cases (e.g. in VCN) it was found advantageous to 
perform the flat fielding partially together with the source fit- 
ting. This was done by adding free parameters to scale inde- 
pendently the measurements made with different detector pix- 
els. This introduces to the fit three additional parameters in the 
case of C200 maps and eight parameters for the C100 maps. 
Determining the flat fielding (multiplicative factor only) locally 
makes it possible to automatically correct for some detector 
drifts. The method is useful when the detector drifts are slow 
compared with the time needed to observe the region surround- 
ing a point source. 

B.2. The fitting procedure 

The point source detection was performed in two steps using 
the surface brightness values and their errors, one value for 
each detector pixel in each raster position. In the first step mea- 
surements more than 0.7<7b g above the local background level 
were flagged as point source candidates. The background level 
and its dispersion, <7bg, were estimated from other measure- 
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merits within a radius which was typically three times the size 
of the detector pixel. 

In the second step a model of a point source and a back- 
ground was fitted to each region surrounding the candidate po- 
sitions. The background was assumed to be constant, since in 
most cases the gradient of the background was small. The free 
parameters of the fit were the source flux density, the two co- 
ordinates of the source position, and the background surface 
brightness. 

The contribution of a point source to the measured flux den- 
sity at different map positions was computed using the foot- 
print matrices of PIA. Footprint matrices provide the fraction 
of the flux detected by each detector pixel. Some 70% of the 
flux from a point source located at the centre of a pixel will 
be detected by this one pixel. The detector beams are approx- 
imately gaussian in their central part but have more extended 
wings. The neighboring pixels will therefore receive slightly 
more flux than predicted by the gaussian approximation, about 
7% or less for source distances exceeding one pixel step (46.0" 
for C100 and 92" for C200). 

The radius of the region used in the fitting procedure was 
typically ~2.2' for maps observed with the detector CI 00 and 
~3.7' for maps done with C200. These radii are large enough 
to make the contribution from the point source small at the edge 
of the fitting region. 

B.3. Completeness and false detections in simulations 

The completeness of the point source detection and the prob- 
ability of false detections were studied with simulated maps 
with raster step equal to the detector size i.e. without any redun- 
dancy. We simulated the dependence of faint sources detection 
on the background noise level. The criterion given above was 
used for selecting the candidate pixels with point source con- 
tribution, the probability limit, P >99%, was used to discard 
uncertain detections. Different values of the ratio 



Sit 

CT bg 



(B.l) 



between source flux density, S v , given in Jy and the back- 
ground rms noise, CTbg, given in units of Jy were also consid- 
ered. The detection rate is about 90% for p — 7.8 while for 
values close to p —5.2 the detection rate drops below 50%. 
The number of false detections as the function of the proba- 
bility limit P was also studied. The number of false detections 
is, as expected, linear with respect to P. However, the number 
of false sources was found to exceed the expected number of 



1-P false detections per map pixel (see Fig. B.la). This was 



taken into account when determining the parameters of source 
detection. One would obviously like to have criteria where the 
number of false detections is equal to the unknown number of 
undetected real sources. 

The flux densities of false detections were found to lie be- 
tween S v w 3.3abg and below S v « 7.0(7b g . The lower limit is 
set partly by the initial selection of candidate pixels that are 
0.7CTb g above the background noise. The spatial distribution 



of the false detections is shown in Fig. B.lb. They are con- 
centrated close to the edges and the corners of the pixels. The 
source density varies roughly as r 2 where r is the distance from 
the pixel centre. This does not, however, help in eliminating 
false detections since the spatial distribution is similar for very 
faint true sources. The bias is present only close to the detection 
limit. 

B.4. Influence of imperfect flat fielding and de-glitching 

The simulations described above, that were used in determin- 
ing suitable detection limits, do not take into account some ef- 
fects that may affect the source counts. For example, imperfect 
flat fielding or de-glitching will lead to higher surface bright- 
ness values in some detector pixels and may increase the prob- 
ability of classifying some random variations as point sources. 
The problem is potentially more serious for the CI 00 detec- 
tor where the deviation of one detector pixel has less influence 
on the overall noise of the surface brightness estimated for the 
region. In those cases where the local flat fielding coefficients 
were estimated together with the source parameters the quality 
of the flat fielding is less important. In order to estimate its ef- 
fect in other cases we simulated CI 00 measurements consisting 
first of normally distributed noise and then scaled upwards the 
values of one detector pixel. The number of false detections did 
not, however, depend strongly on the imperfections of the flat 
fielding and was in fact smaller for the tested range where one 
detector pixel deviated less than 2a. If more than one pixel de- 
viates from the mean the effect is further reduced. We conclude 
that imperfect flat fielding will not cause significant errors in 
the source counts. 



Appendix C: Methods in the examination of SRD data 

The main procedures used in the source extraction were based 
on AAP (Astrophysical Applications Data) data products re- 
duced with PIA. The corresponding SRD (Signal per Ramp 
Data) files were used for visual inspection of the data. Glitches 
are visible at the SRD level as a sudden high signal value fol- 
lowed by gradually decreasing tail or, in the case of several 
glitches close in time, as unusually large noise. The number of 
ramps per raster position varies in our data from 6 to over 40. 
In the case of the CI 00 detector, which is more affected by the 
glitches, the number of ramps was always sufficient to see the 
characteristic features of the glitches. 

The SRD data close to the potential sources were inspected 
visually in order to check whether the detections were possibly 
due to glitches and not to real point sources. The raster posi- 
tion and the detector pixel closest to the fitted source position 
was determined and with the selected raster position in the cen- 
tre the SRD data for five consecutive raster position and for all 
detector pixels were plotted (see Fig. ^|). Both the variations in 
the signal level between raster positions and between the differ- 
ent detector pixels were examined by eye and the measurement 
closest to the source position was checked for signs of glitches. 
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Fig. B.l. a The number of false detections per map pixel in 
simulated measurements as a function of the expected number 
1-P. The confidence level P is used to discard uncertain de- 
tections. At high confidence levels the actual number of false 
detections exceeds the expected number (dashed line), b Posi- 
tions of false detections inside a map pixel for non-overlapping 
rasters. Similar distribution exists for very weak sources 



A quality flag between and 4 was given for the source and 
these are shown in Table 5.1 as flag "q". The intended scale is: 



- 0: the median signal is clearly affected by glitches and/or 
the noise is clearly too high for reliable determination of 
the median signal 

- 1 : there are clear glitches that may have affected the median 
signal and/or large noise makes the determination of the 
median uncertain 

- 2: no clear signs of a source, no significant glitches 

- 3: a possible source 

- 4: a clear detection 

The difference in the signal levels between different raster 
positions and between different detector pixels were also used 
in the scrutinization. Therefore, even if the signals were seri- 
ously affected by glitches the source could be classified in e.g. 
class 3 provided that the true signal level could be estimated to 
be high enough. 

The classification is, of course, subjective. It is also strongly 
biased towards sources that happen to lie in the centre of some 
detector pixel and are therefore visible only in one measure- 
ment. Furthermore, if the source is situated between two raster 
lines the previously described SRD vs. time plots do not show 
all data relevant for the classification. 

The main benefit from the eyeball inspection is that we can 
recognize potentially false detections (classes and 1) that are 
due to detector glitches. There are, however, only a couple of 
such possibly false detections in our source sample. This shows 
that most false detections were avoided in our source detection 
procedure. 

The source detections were tested also by applying to the 
SRD data methods reminisc ent of the procedures used in the 
ELAIS project (Surace et al. |1999| ). 

For each source we determined the raster position and the 
detector pixel, po closest to the potential source. Using the sig- 
nal values sf from all detector pixels p from five consecutive 
raster positions (i=-2. . . 2) centered on the selected position we 
calculated the following values: 



1. 



s o° / s o ' P 7^ Po subtracted by the median of the correspond- 
ing values at the five raster positions 

2. same as (1.) but instead of the median we use a prediction 
from 2nd degree fit to the five raster positions 

3. so subtracted by the median of sf°, i=-2. . . 2 

4. same as previous but using prediction from a 2nd degree fit 
instead of the median 

5. Sq° divided by the median of other pixels in the current 
raster position subtracted by the median of the correspond- 
ing values in all five raster positions 

6. same as previous but using a 2nd degree fit instead of the 
median of five raster positions 

In the case of CI 00 we calculated therefore altogether 20 
values and in the case of C200 10 values. The number of val- 
ues above 1.5(7 level was counted and is shown in Table 5.1 as 



flag "c". It must be emphasized that the method is applied only 
to the SRD data closest to the previously determined source 
positions and it is therefore most sensitive to sources that were 
close to the centre of some detector pixel. Thus a low value of 
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"c" does not disqualify a source detection obtained with our 
detection procedure. 
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Table 4: IRAS PSC and FSC sources located within the ISOPHOT maps and the 
data on our FIR detections. First columns indicates the IRAS source catalog, 
columns (2)-(10) give the positions and flux densities of the IRAS sources as well 
as the sizes of the error ellipses. For ISOPHOT observations we quote for non- 
tections upper limits for the flux densities (^2.5er limit). For detected sources 
give at each wavelength separately the distance from the IRAS position 
pLumns (11), (13), (15) and (17)) 
IRAS 
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(2000.0) 
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(") 
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68 


0.34 


0.84 
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8 
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16 


1 


68 
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<1.00 


33 


6 


68 
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26 


3 


68 


0.75 
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59 


13 


87 
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<0.99 


27 


7 
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0.23 
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30 


9 
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6 


2 
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28 


11 
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Table 5: Galaxies and quasars in the mapped regions from the Simbad database. 
The list contains also other sources when they are located close to our FIR 
detections. The columns are: (1) the name of the source, (2)-(3) coordinates, (4) 
the classification of the source, (5)-(12) data based on ISOPHOT observations: 
positional distances from Simbad source and our estimates for the FIR flux 
densities. In case no FIR sources were detected close to the Simbad positions 
estimated upper limits for the flux densities are given 
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1 Source close to the border of the mapped area 

2 Galaxy Z 1339.5+4108 does not exist in our list of source detections. It lies 
at the edge of the large 180yum map i.e. it has been observed at only one 
wavelength. We did include into the source list those sources that could be 
detected in maps obtained by dividing the data of this 180^m map at the SRD 
level into two halves. At the SRD level each measurement of a map position 
contained just 6 signal values i.e. after the division each of the two maps had 

3 signals per position. At the position of the galaxy five of the six signals were 
high but the last one was below the background level. The galaxy was therefore 
clearly detected in the first map but the confidence obtained from the second 
remained below our detection threshold. The quoted flux density estimate is 
based on all of the 180/im data. 
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